Predicting enteric methane emission in lactating Holsteins based on reference methane data collected by the GreenFeed system.

Methane emission is not included in the current breeding goals for dairy cattle mainly due to the expense and difficulty in obtaining sufficient data to generate accurate estimates of the relevant traits. While several models have been developed to predict methane emission from milk spectra using reference methane data obtained by the respiration chamber, SF6 and sniffer methods, the prediction of methane emission from milk mid-infrared (MIR) spectra using reference methane data collected by the GreenFeed system has not yet been explored. Methane emission was monitored for 151 cows using the GreenFeed system. Prediction models were developed for daily and average (for the trial period of 12 or 14 days) methane production (g/d), yield (g/kg DM intake (DMI)) and intensity (g/kg of fat- and protein-corrected milk) using partial least squares regression. The predictions were evaluated in 100 repeated validation cycles, where animals were randomly partitioned into training (80%) and testing (20%) populations for each cycle. The best performing model was observed for average methane intensity using MIR, parity and DMI with validation coefficient of determination (R2val) and RMSE of prediction of 0.66 and 4.7 g/kg of fat- and protein-corrected milk, respectively. The accuracy of the best models for average methane production and average methane yield were poor (R2val = 0.28 and 0.12, respectively). A lower accuracy of prediction was observed for methane intensity and production (R2val = 0.42 and 0.17) when daily records were used while prediction for methane yield was comparable to that for average methane yield (R2val = 0.16). Our results suggest the potential to predict methane intensity with moderate accuracy. In this case, prediction models for average methane values were generally better than for daily measures when using the GreenFeed system to obtain reference methane emission measurements.

Breeding for reduced methane emission and feed-efficient Holstein cows: An international response.

Selecting for lower methane (CH4) emitting animals is one of the best approaches to reduce CH4 given that genetic progress is permanent and cumulative over generations. As genetic selection requires a large number of animals with records and few countries actively record CH4, combining data from different countries could help to expedite accurate genetic parameters for CH4 traits and build a future genomic reference population. Additionally, if we want to include CH4 in the breeding goal, it is important to know the genetic correlations of CH4 traits with other economically important traits. Therefore, the aim of this study was first to estimate genetic parameters of 7 suggested methane traits, as well as genetic correlations between methane traits and production, maintenance, and efficiency traits using a multicountry database. The second aim was to estimate genetic correlations within parities and stages of lactation for CH4. The third aim was to evaluate the expected response of economically important traits by including CH4 traits in the breeding goal. A total of 15,320 methane production (MeP, g/d) records from 2,990 cows belonging to 4 countries (Canada, Australia, Switzerland, and Denmark) were analyzed. Records on dry matter intake (DMI), body weight (BW), body condition score, and milk yield (MY) were also available. Additional traits such as methane yield (MeY; g/kg DMI), methane intensity (MeI; g/kg energy-corrected milk), a genetic standardized methane production, and 3 definitions of residual methane production (g/d), residual feed intake, metabolic BW (MBW), BW change, and energy-corrected milk were calculated. The estimated heritability of MeP was 0.21, whereas heritability estimates for MeY and MeI were 0.30 and 0.38, and for the residual methane traits heritability ranged from 0.13 to 0.16. Genetic correlations between different methane traits were moderate to high (0.41 to 0.97). Genetic correlations between MeP and economically important traits ranged from 0.29 (MY) to 0.65 (BW and MBW), being 0.41 for DMI. Selection index calculations showed that residual methane had the most potential for inclusion in the breeding goal when compared with MeP, MeY, and MeI, as residual methane allows for selection of low methane emitting animals without compromising other economically important traits. Inclusion of residual feed intake in the breeding goal could further reduce methane, as the correlation with residual methane is moderate and elicits a favorable correlated response. Adding a negative economic value for methane could facilitate a substantial reduction in methane emissions while maintaining an increase in milk production.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing.

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled.

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

Development and validation of a small SNP panel for feed efficiency in beef cattle.

The objective of this study was to develop and validate a customized cost-effective single nucleotide polymorphism (SNP) panel for genetic improvement of feed efficiency in beef cattle. The SNPs identified in previous association studies and through extensive analysis of candidate genomic regions and genes, were screened for their functional impact and allele frequency in Angus and Hereford breeds used as validation candidates for the panel. Association analyses were performed on genotypes of 159 SNPs from new samples of Angus (n = 160), Hereford (n = 329), and Angus-Hereford crossbred (n = 382) cattle using allele substitution and genotypic models in ASReml. Genomic heritabilities were estimated for feed efficiency traits using the full set of SNPs, SNPs associated with at least one of the traits (at P ≤ 0.05 and P < 0.10), as well as the Illumina bovine 50K representing a widely used commercial genotyping panel. A total of 63 SNPs within 43 genes showed association (P ≤ 0.05) with at least one trait. The minor alleles of SNPs located in the GHR and CAST genes were associated with decreasing effects on residual feed intake (RFI) and/or RFI adjusted for backfat (RFIf), whereas minor alleles of SNPs within MKI67 gene were associated with increasing effects on RFI and RFIf. Additionally, the minor allele of rs137400016 SNP within CNTFR was associated with increasing average daily gain (ADG). The SNPs genotypes within UMPS, SMARCAL, CCSER1, and LMCD1 genes showed significant over-dominance effects whereas other SNPs located in SMARCAL1, ANXA2, CACNA1G, and PHYHIPL genes showed additive effects on RFI and RFIf. Gene enrichment analysis indicated that gland development, as well as ion and cation transport are important physiological mechanisms contributing to variation in feed efficiency traits. The study revealed the effect of the Jak-STAT signaling pathway on feed efficiency through the CNTFR, OSMR, and GHR genes. Genomic heritability using the 63 significant (P ≤ 0.05) SNPs was 0.09, 0.09, 0.13, 0.05, 0.05, and 0.07 for ADG, dry matter intake, midpoint metabolic weight, RFI, RFIf, and backfat, respectively. These SNPs contributed to genetic variation in the studied traits and thus can potentially be used or tested to generate cost-effective molecular breeding values for feed efficiency in beef cattle.

Comparative transcriptomic analysis of porcine peripheral blood reveals differentially expressed genes from the cytokine-cytokine receptor interaction pathway related to health status.

While some research has looked into the host genetic response in pigs challenged with specific viruses or bacteria, few studies have explored the expression changes of transcripts in the peripheral blood of sick pigs that may be infected with multiple pathogens on farms. In this study, the architecture of the peripheral blood transcriptome of 64 Duroc sired commercial pigs, including 18 healthy animals at entry to a growing facility (set as a control) and 23 pairs of samples from healthy and sick pen mates, was generated using RNA-Seq technology. In total, 246 differentially expressed genes were identified to be specific to the sick animals. Functional enrichment analysis for those genes revealed that the over-represented gene ontology terms for the biological processes category were exclusively immune activity related. The cytokine-cytokine receptor interaction pathway was significantly enriched. Nine functional genes from this pathway encoding members (as well as their receptors) of the interleukins, chemokines, tumor necrosis factors, colony stimulating factors, activins, and interferons exhibited significant transcriptional alteration in sick animals. Our results suggest a subset of novel marker genes that may be useful candidate genes in the evaluation and prediction of health status in pigs under commercial production conditions.

Genome-wide association study for lactation persistency, female fertility, longevity, and lifetime profit index traits in Holstein dairy cattle.

Female fertility in Holstein cattle can decline when intense genetic selection is placed on milk production. One approach to improving fertility is to identify the genomic regions and variants affecting fertility traits and then incorporate this knowledge into selection decisions. The objectives of this study were to identify or refine the positions of the genomic regions associated with lactation persistency, female fertility traits (age at first service, cow first service to conception, heifer and cow nonreturn rates), longevity traits (herd life, indirect herd life, and direct herd life), and lifetime profit index in the North American Holstein dairy cattle population. A genome-wide association study was performed for each trait, using a single SNP (single nucleotide polymorphism) regression mixed linear model and imputed high-density panel (777k) genotypes. No associations were identified for fertility traits. Several peak regions were detected for lifetime profit index, lactation persistency, and longevity. The results overlap with previous findings and identify some novel regions for lactation persistency. Previously proposed causative and candidate genes supported by this work include DGAT1, GRINA, and CPSF1, whereas new candidate genes are SLC2A4RG and THRB. Thus, the chromosomal regions identified in this study not only confirm several previous findings but also highlight new regions that may contribute to genetic variation in lactation persistency and longevity-associated traits in dairy cattle.

Accuracy of genomic predictions for feed efficiency traits of beef cattle using 50K and imputed HD genotypes.

The accuracy of genomic predictions can be used to assess the utility of dense marker genotypes for genetic improvement of beef efficiency traits. This study was designed to test the impact of genomic distance between training and validation populations, training population size, statistical methods, and density of genetic markers on prediction accuracy for feed efficiency traits in multibreed and crossbred beef cattle. A total of 6,794 beef cattle data collated from various projects and research herds across Canada were used. Illumina BovineSNP50 (50K) and imputed Axiom Genome-Wide BOS 1 Array (HD) genotypes were available for all animals. The traits studied were DMI, ADG, and residual feed intake (RFI). Four validation groups of 150 animals each, including Angus (AN), Charolais (CH), Angus-Hereford crosses (ANHH), and a Charolais-based composite (TX) were created by considering the genomic distance between pairs of individuals in the validation groups. Each validation group had 7 corresponding training groups of increasing sizes ( = 1,000, 1,999, 2,999, 3,999, 4,999, 5,998, and 6,644), which also represent increasing average genomic distance between pairs of individuals in the training and validations groups. Prediction of genomic estimated breeding values (GEBV) was performed using genomic best linear unbiased prediction (GBLUP) and Bayesian method C (BayesC). The accuracy of genomic predictions was defined as the Pearson's correlation between adjusted phenotype and GEBV (), unless otherwise stated. Using 50K genotypes, the highest average achieved in purebreds (AN, CH) was 0.41 for DMI, 0.34 for ADG, and 0.35 for RFI, whereas in crossbreds (ANHH, TX) it was 0.38 for DMI, 0.21 for ADG, and 0.25 for RFI. Similarly, when imputed HD genotypes were applied in purebreds (AN, CH), the highest average was 0.14 for DMI, 0.15 for ADG, and 0.14 for RFI, whereas in crossbreds (ANHH, TX) it was 0.38 for DMI, 0.22 for ADG, and 0.24 for RFI. The of GBLUP predictions were greatly reduced with increasing genomic average distance compared to those from BayesC predictions. The results indicate that 50K genotypes, used with BayesC, are more effective for predicting GEBV in purebred cattle. Imputed HD genotypes found utility when dealing with composites and crossbreds. Formulation of a fairly large training set for genomic predictions in beef cattle should consider the genomic distance between the training and target populations.

Transcriptomic analysis by RNA sequencing reveals that hepatic interferon-induced genes may be associated with feed efficiency in beef heifers.

In beef cattle, production feedstuffs are the largest variable input cost. Beef cattle also have a large carbon footprint, raising concern about their environmental impact. Unfortunately, only a small proportion of dietary energy is directed toward protein deposition and muscle growth whereas the majority supports body maintenance. Improving feed efficiency would, therefore, have important consequences on productivity, profitability, and sustainability of the beef industry. Various measures of feed efficiency have been proposed to improve feed utilization, and currently, residual feed intake (RFI) is gaining popularity. However, the cost associated with measuring RFI and the limited knowledge of the biology underlying improved feed efficiency make its adoption prohibitive. Identifying molecular mechanisms explaining divergence in RFI in beef cattle would lead to the development of early detection methods for the selection of more efficient breeding stock. The objective of this study was to identify hepatic markers of metabolic feed efficiency in replacement beef heifers. A group of 87 heifers were tested for RFI adjusted for off-test backfat thickness (RFIfat). Preprandial liver biopsies were collected from 10 high- and 10 low-RFIfat heifers (7 Hereford­Aberdeen Angus and 3 Charolais­Red Angus­Main Anjou per group) and gene expression analysis was performed using RNA sequencing and quantitative real-time PCR. The heifers used in this study differed in RFIfat averaging 0.438 vs. ­0.584 kg DM/d in high- and low-RFIfat groups, respectively. As expected, DMI was correlated with RFIfat and ADG did not differ between high- and low-RFIfat heifers. Through a combination of whole transcriptome and candidate gene analyses, we identified differentially expressed genes involved in inflammatory processes including hemoglobin ß (HBB), myxovirus resistance 1 interferon-inducible protein p78 (MX1), ISG15 ubiquitin-like modifier (ISG15), hect domain and RLD 6 (HERC6), and interferon-induced protein 44 (IFI44) whose mRNA abundance was lower (HBB) or higher (MX1, ISG15, HERC6, and IFI44) in low-RFIfat heifers. These genes have been shown to be directly or indirectly modulated by interferon signaling and involved with innate immunity. Our results suggest that more efficient heifers respond differently to hepatic proinflammatory stimulus, potentially expending less energy toward combating systemic inflammation and redirecting nutrients toward growth and protein accretion.

Candidate genes and single nucleotide polymorphisms associated with variation in residual feed intake in beef cattle.

The candidate gene approach was used to identify genes associated with residual feed intake (RFI) in beef steers. The approach uses prior knowledge of gene functions to predict their biological role in the variation observed in a trait. It is suited to identify genes associated with complex traits where each gene has a relatively small effect. First, positional candidate genes were identified within the genomic positions of previously reported QTL associated with component traits related to RFI such as dry matter intake (DMI), growth, feed conversion ratio (FCR), average daily gain (ADG), and energy balance. Secondly, the positional candidate genes were prioritized into functional candidate genes according to their biological functions and their relationship with the biological processes associated with RFI including carbohydrate, fat and protein metabolism, thermoregulation, immunity and muscle activity. Single nucleotide polymorphisms (SNPs) located within the functional candidate genes were identified using mRNA sequences and prioritized into functional classes such as non-synonymous (nsSNP), synonymous (sSNP) or intronic SNP. A total of 117 nsSNP were considered as functional SNP and genotyped in steers at the University of Alberta ranch in Kinsella. Multiple marker association analysis in ASReml was performed using RFI data obtained from 531 beef steers. Twenty-five SNP were significantly associated with RFI (P < 0.05) accounting for 19.7% of the phenotypic variation. Using SIFT program to predict the effect of the SNP on the function of the corresponding protein, 3 of the 25 SNP were predicted to cause a significant effect on protein function (P < 0.05). One of the 3 SNP was located in the GHR gene and was also associated with a significant effect on the tertiary structure of the GHR protein (P < 0.05) as modeled using SWISSModel software. Least square means for each genotype were estimated and an over-dominance effect was observed for the SNP located in the GHR, CAST, ACAD11 and UGT3A1 genes. In addition, 2 other SNP showed a dominance effect and 3 genes had an additive effect. Gene network analysis performed in Ingenuity pathway analysis (IPA) software (Ingenuity Systems, www.ingenuity.com) indicated that the significant genes were involved in biological pathways such as lipid, protein and energy metabolism, electron transport and membrane signaling. The genes in this study, if validated in other beef cattle populations, may be useful for marker assisted selection for feed efficiency.

Whole-genome QTL scan for ultrasound and carcass merit traits in beef cattle using Bayesian shrinkage method.

Fine mapping of quantitative trait loci (QTL) for 16 ultrasound measurements and carcass merit traits that were collected from 418 hybrid steers was conducted using 1207 SNP markers covering the entire genome. These SNP markers were evaluated using a Bayesian shrinkage estimation method and the empirical critical significant thresholds (α = 0.05 and α = 0.01) were determined by permutation based on 3500 permuted datasets for each trait to control the genome-wide type I error rates. The analyses identified a total of 105 QTLs (p < 0.05) for seven ultrasound measure traits including ultrasound backfat, ultrasound marbling and ultrasound ribeye area and 113 QTLs for seven carcass merit traits of carcass weight, grade fat, average backfat, ribeye area, lean meat yield, marbling and yield grade. Proportion of phenotypic variance accounted for by a single QTL ranged from 0.06% for mean ultrasound backfat to 4.83% for carcass marbling (CMAR) score, while proportion of the phenotypic variance accounted for by all significant (p < 0.05) QTL identified for a single trait ranged from 4.54% for carcass weight to 23.87% for CMAR.

Identification of candidate markers on bovine chromosome 14 (BTA14) under milk production trait quantitative trait loci in Holstein.

The objective of this study was to identify single-nucleotide polymorphisms using a bovine chromosome 14 high-density SNP panel after accounting for the effect of DGAT1. Linkage disequilibrium information and sire heterozygosity were used to select markers for linkage analysis on bovine chromosome 14 for milk production traits in 321 Holstein animals. Results show putative milk peaks at 42 and 61 cM, both at p<0.10, a fat yield peak at 42 and 63 cM, both at p<0.05; a protein yield peak at 42 (p<0.01) and 84 cM (p<0.05); fat per cent peaks at 3 (p<0.01) and 29 cM (p<0.05), and a protein per cent peak at 4 cM (p<0.05). Once quantitative trait loci positions were established, allele substitution effects for all markers were evaluated using the same statistical model. Overlaying information between quantitative trait loci (QTL) and allele effect analysis enabled the identification (p<0.01) of 20 SNPs under the milk yield QTL, 2 under both of the fat yield peaks, 8 and 9 under the protein yield peaks, 2 and 6 for the fat per cent peaks and 5 for the protein per cent peak. One SNP in particular, ss61514555:A>C, showed association with 3 of the 5 traits: milk (p=1.59E-04), fat (p=6.88E-05) and protein yields (p=5.76E-05). Overall, combining information from linkage disequilibrium, sire heterozygosity and genetic knowledge of traits enabled the characterization of additional markers with significant associations with milk production traits.

Accuracy of genomic breeding values for residual feed intake in crossbred beef cattle.

The benefit of using genomic breeding values (GEBV) in predicting ADG, DMI, and residual feed intake for an admixed population was investigated. Phenotypic data consisting of individual daily feed intake measurements for 721 beef cattle steers tested over 5 yr was available for analysis. The animals used were an admixed population of spring-born steers, progeny of a cross between 3 sire breeds and a composite dam line. Training and validation data sets were defined by randomly splitting the data into training and testing data sets based on sire family so that there was no overlap of sires in the 2 sets. The random split was replicated to obtain 5 separate data sets. Two methods (BayesB and random regression BLUP) were used to estimate marker effects and to define marker panels and ultimately the GEBV. The accuracy of prediction (the correlation between the phenotypes and GEBV) was compared between SNP panels. Accuracy for all traits was low, ranging from 0.223 to 0.479 for marker panels with 200 SNP, and 0.114 to 0.246 for marker panels with 37,959 SNP, depending on the genomic selection method used. This was less than accuracies observed for polygenic EBV accuracies, which ranged from 0.504 to 0.602. The results obtained from this study demonstrate that the utility of genetic markers for genomic prediction of residual feed intake in beef cattle may be suboptimal. Differences in accuracy were observed between sire breeds when the random regression BLUP method was used, which may imply that the correlations obtained by this method were confounded by the ability of the selected SNP to trace breed differences. This may also suggest that prediction equations derived from such an admixed population may be useful only in populations of similar composition. Given the sample size used in this study, there is a need for increased feed intake testing if substantially greater accuracies are to be achieved.

A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls.

The detection and mapping of genetic markers linked to quantitative trait loci (QTL) can be utilized to enhance genetic improvement of livestock populations. With the completion of the bovine genome sequence assembly, single nucleotide polymorphisms (SNP) assays spanning the whole bovine genome and research work on large scale identification, validation and analysis of genotypic variation in cattle has become possible. The objective of the present study was to perform a whole genome scan to identify and map QTL affecting milk production traits and somatic cell scores using linkage disequilibrium (LD) regression and 1536 SNP markers. Three and 18 SNP were found to be associated with only milk yield (MY) at a genome and chromosome wise significance (p < 0.05) level respectively. Among the 21 significant SNP, 16 were in a region reported to have QTL for MY in other dairy cattle populations and while the rest five were new QTL finding. Four SNP out of 21 are significant for the milk production traits (MY, fat yield, protein yield (PY), and milk contents) in the present study. Six and nine SNP were associated with PY at a genome and chromosome wise significant (p < 0.05) level respectively. Three and 17 SNP were found to be associated with FY at a genome and chromosome wise significant (p < 0.05) level. Five and seven SNP were mapped with somatic cell score at a genome and chromosome wise significant (p < 0.05) level respectively. The results of this study have revealed QTL for MY, PY, protein percentage, FY, fat percentage, somatic cell score and persistency of milk in the Canadian dairy cattle population. The chromosome regions identified in this study should be further investigated to potentially identify the causative mutations underlying the QTL.

Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle.

The objective of this study was to quantify the extent of linkage disequilibrium (LD) on bovine chromosomes 19 and 29 and to study the pattern of selection signatures in beef and dairy breeds (Angus and Holstein) of Bos taurus. The extent of LD was estimated for 370 and 186 single nucleotide polymorphism markers on BTA19 and 29 respectively using the square of the correlation coefficient (r(2)) among alleles at pairs of loci. A comparison of the extent of LD found that the decline of LD followed a similar pattern in both breeds. We observed long-range LD and found that LD dissipates to background levels at a locus separation of about 20 Mb on both chromosomes. Along each chromosome, patterns of LD were variable in both breeds. We find that a minimum of 30 000 informative and evenly spaced markers would be required for whole-genome association studies in cattle. In addition, we have identified chromosomal regions that show some evidence of selection for economically important traits in Angus and Holstein cattle. The results of this study are of importance for the design and application of association studies.

A whole-genome scan to map quantitative trait loci for conformation and functional traits in Canadian Holstein bulls.

Genetic improvement of livestock populations can be achieved through detection and mapping of genetic markers linked to quantitative trait loci (QTL). With the completion of the bovine genome sequence assembly, single nucleotide polymorphism (SNP) assays spanning the whole bovine genome and research work on large-scale identification, validation, and analysis of genotypic variation in cattle has become possible. A total of 462 Canadian Holstein Bulls were used to test the association between SNP and QTL. Single locus linkage disequilibrium regression model was implemented to perform a whole genome scan to identify and map QTL affecting conformation and functional traits. One thousand five hundred thirty-six SNP markers from introns and exons of potential QTL regions for economically important traits across the bovine genome were selected for association analysis. A total of 45 and 151 SNP were found to be associated with 17 conformation and functional traits at a genome- and chromosome-wise significance level, respectively. Among the 196 significant SNP, 169 of them are newly detected in this study, whereas 27 of them have been reported in previous literature and 161 of these were located in genes and are worth further investigating to potentially identify the causative mutations underlying the QTL. The single locus linkage disequilibrium regression method using SNP marker genotypes has proven to be a successful methodology for detecting and mapping QTL in dairy cattle populations.

Construction of bovine whole-genome radiation hybrid and linkage maps using high-throughput genotyping.

High-density whole-genome maps are essential for ordering genes or markers and aid in the assembly of genome sequence. To increase the density of markers on the bovine radiation hybrid map, and hence contribute to the assembly of the bovine genome sequence, an Illumina BeadStation was used to simultaneously type large numbers of markers on the Roslin-Cambridge 3000 rad bovine-hamster whole-genome radiation hybrid panel (WGRH3000). In five multiplex reactions, 6738 sequence tagged site (STS) markers were successfully typed on the WGRH3000 panel DNA. These STSs harboured SNPs that were developed as a result of the bovine genome sequencing initiative. Typically, the most time consuming and expensive part of creating high-density radiation hybrid (RH) maps is genotyping the markers on the RH panel with conventional approaches. Using the method described in this article, we have developed a high-density whole-genome RH map with 4690 loci and a linkage map with 2701 loci, with direct comparison to the bovine whole-genome sequence assembly (Btau_2.0) in a fraction of the time it would have taken with conventional typing and genotyping methods.

COMBOSA3D: combining sequence alignments with three-dimensional structures.

UNLABELLED: COMBOSA3D is a program that allows sequence conservation to be viewed in its proper three-dimensional environment. It superimposes sequence alignment information onto a protein structure using a customizable color scheme, which is also applied to a textual sequence alignment for reference. AVAILABILITY: The program can be tested at http://www.bioinformatics.org/combosa3d/, and the source code is freely available.

PepTool and GeneTool: platform-independent tools for biological sequence analysis.

Although we are unable to discuss all of the functionality available in PepTool and GeneTool, it should be evident from this brief review that both packages offer a great deal in terms of functionality and ease-of-use. Furthermore, a number of useful innovations including platform-independent GUI design, networked parallelism, direct internet connectivity, database compression, and a variety of enhanced or improved algorithms should make these two programs particularly useful in the rapidly changing world of biological sequence analysis. More complete descriptions of the programs, algorithms and operation of PepTool and GeneTool are available on the BioTools web site (www.biotools.com), in the associated program user manuals and in the on-line Help pages.